clc; clear all;
H= 19e6;
Rz=6400e3;
Tobr=12*60*60;
Tmod=(Tobr/2);%vremia modelirovania
wobr=2*pi/Tobr;
c_light=3e8;
f_0=1602e6;%chastota
dt=5*50;
t=0:dt:Tmod;

%vector potrebitelia
x=[0;Rz]; %koordinata potribitelia



for i=1:length(t)
rsv=(H+Rz)*...
    [cos(wobr*t(i));
    sin(wobr*t(i))];
plot( [0 rsv(1)],[0 rsv(2)]);
xlim([-H-Rz;H+Rz])
ylim([-Rz;H+Rz])
axis equal
hold on
pause (0.1);
end







